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Abstract 



D0 has measured the inclusive production cross section of W and Z bosons 
in a sample of 13 pb _1 of data collected at the Fermilab Tevatron. The cross 
sections, multiplied by their leptonic branching fractions, for production in pp 
collisions at y/s = 1.8 TeV are a w B(W -► ev) = 2.36 ± 0.02 ± 0.08 ± 0.13 nb, 
a w B(W -► nv) = 2.09 ± 0.06 ± 0.22 ± 0.11 nb, a z B(Z -► e + e~) = 0.218 ± 
0.008 ± 0.008 ± 0.012 nb, and a z B(Z -► n + n~) = 0.178 ± 0.022 ± 0.021 ± 
0.009 nb, where the first uncertainty is statistical and the second systematic; 
the third reflects the uncertainty in the integrated luminosity. For the com- 
bined electron and muon analyses, we find <7jyB(VF — ► lu)/azB(Z — > l + l~) = 
10.90 ± 0.52. Assuming standard model couplings, we use this result to de- 
termine the width of the W boson, and obtain T(W) = 2.044 ± 0.097 GeV. 
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I. INTRODUCTION 



Measurement of the production cross sections multiplied by the leptonic branching frac- 
tions (crB) for W and Z bosons can be used to test predictions of QCD for W and Z 
production, and to extract the width of the W boson. Previous measurements of these cross 
sections have been made at y/s = 630 GeV by the UA1 (lj and UA2 experiments and at 
y/s = 1.8 TeV by CDF 0-0. The results reported in this paper are from the D0 detector, 
operating at = 1.8 TeV, and have been summarized previously in Ref. ||. 

Precise determination of the total widths of the W and Z bosons provides an important 
test of the standard model because these widths are sensitive to new (and possibly unde- 
tected) decay modes. The total width of the Z boson is known to a precision of 0.3% || 
which places strong constraints on the existence of any new particles that can contribute 
to decays to neutrals. Our knowledge of the total width of the W boson is an order of 
magnitude less precise, and the corresponding limits on charged weak decays are much less 
stringent. It is therefore important to improve the measurement of the width of the W 
boson as a means of searching for any unexpected PF-boson decay modes. 

We determine the width of the W boson indirectly by using the ratio of the measured 
W and Z boson aB values 

a w B(W -> lu) 
azB(Z^ll) ' 

where / corresponds to e or /i, aw and az are the inclusive cross sections for W and Z boson 
production, a(pp — > W + X) and a(pp — > Z + X), and B(W — > lu) and B{Z — > 11) are the 
leptonic branching fractions of the W and Z bosons. We extract B{W — > lu) from the above 
ratio using a theoretical prediction for aw/cz, and the precise measurement of B(Z — > 11) 
from LEP. We then combine B(W — > lu) with the leptonic partial width T(W — > lu) to 
obtain the total width of the W boson, r(W). 

Many of the systematic uncertainties (both experimental and theoretical) that affect the 
determination of the individual cross sections ayi/P^VF — > lu) and azBiZ — > 11) cancel when 
calculating the ratio R. At the present time, R gives the best determination of T(W / ); direct 
measurements from fits to the tail of the transverse mass distribution of the W boson are 



currently four times less precise |TTJ, but require fewer standard model assumptions. 

In this paper, we report results of the measurement of the W and Z production cross 
sections, and the extraction of r(W^), using data collected in the first collider run of the D0 
detector starting in August 1992 and ending in June 1993. During the run, the Tevatron 
operated with a typical instantaneous luminosity of 4.0 x 10 30 cm _2 s _1 and a peak luminosity 
of 9.7 x 10 30 cm _2 s _1 . D0 recorded to tape a total of ~ 13 pb _1 of data. 



II. THE D0 DETECTOR 

D0 is a multipurpose detector designed to study pp collisions at the Fermilab Tevatron 
Collider. It consists of three primary components: a nonmagnetic central tracking system, 
a nearly hermetic uranium/liquid- argon calorimeter, and a magnetic muon spectrometer. A 
cutaway view of the detector is shown in Fig. p]. A full description can be found in Ref. [|i~2" 
below we give details of the detector relevant to this analysis. 
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FIG. 1. Cutaway isometric view of the D0 detector. 

A. Central Tracking System 

The central tracking system consists of four detector subsystems: a vertex drift cham- 
ber (VTX), a transition radiation detector (TRD), a central drift chamber (CDC), and 
two forward drift chambers (FDC). The system provides charged-particle tracking over the 
pseudorapidity region \r)\ < 3.2, where 77 = tanh _1 (cos 9), 9 is the polar angle, and is 
the azimuthal angle. Trajectories of charged particles are measured with a resolution of 
2.5 mrad in and 28 mrad in 9. From these measurements, the position of the interac- 
tion vertex along the beam direction (z) can be determined with a typical resolution of 8 
mm. The central tracking system also measures the ionization of tracks, and can be used to 
distinguish single charged particles and e + e~ pairs from photon conversions. 

B. Calorimeter 

Surrounding the central tracking system is the calorimeter, which is divided into three 
parts: a central calorimeter (CC) and two end calorimeters (EC). They each consist of an 
inner electromagnetic (EM) section, a fine hadronic (FH) section, and a coarse hadronic (CH) 
section, housed in a steel cryostat. The intercryostat detector (ICD) consists of scintillator 
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tiles inserted in the space between the EC and CC cryostats. The ICD improves the energy 
resolution for jets that straddle two cryostats. The calorimeter covers the range \rj\ < 4.2. 

Each EM section is 21 radiation lengths deep, and is divided into four longitudinal 
segments (layers). The hadronic sections are 7-9 nuclear interaction lengths deep, and are 
divided into four (CC) or five (EC) layers. The calorimeter is transversely segmented into 
pseudoprojective towers of Arj x A0 = 0.1 x 0.1. The third layer of the EM calorimeter, 
in which the maximum energy deposition of EM showers is expected, is segmented twice 
as finely into cells with Arj x A<p = 0.05 x 0.05. With this fine segmentation, the position 
resolution for electrons above 50 GeV in energy is about 2.5 mm. The energy resolution is 
a(E)/E = 15%/ J E (GeV) © 0.4% for electrons. For charged pions the resolution is about 



50%/ v /E(GeV), and for jets about 80%/yJE(GeV) pj. From minimum-bias data, for the 
imbalance in transverse momentum, or "missing EV' (see Sec. [Ill C| ), or 1/t, the resolution 
for each component (If x and 1/J y ) is 1.08 GeV + 0.019(T,Et), where T,Et is the scalar sum 
of the transverse energies over all calorimeter cells. 

The readout of the individual calorimeters cells are subject to zero suppression. They 
are read out only if the signal is outside of a two-standard-deviations window centered on 
the mean of the noise. 



C. Muon Spectrometer 

Outside the calorimeter, there are muon detection systems covering \r]\ < 3.3. Since 
muons from W and Z boson decays populate predominantly the central region, this work 
uses only the wide angle muon spectrometer (WAMUS), which consists of four planes of 
proportional drift tubes (PDT) in front of magnetized iron toroids with a magnetic field of 
1.9 T, and two groups of three planes each of proportional drift tubes behind the toroids. 
The magnetic field lines and the wires in the drift tubes are oriented transversely to the 
beam direction. The WAMUS covers the region < 1.7 over the entire azimuth, with the 
exception of the central region below the calorimeter (\r]\ < 1, 225° < <p < 315°), where the 
inner layer is missing to make room for the support structure of the calorimeter. 

The total material in the calorimeter and iron toroids varies between 13 and 19 in- 
teraction lengths, making background in the muon chambers from hadronic punchthrough 
negligible. The D0 detector is significantly more compact than previous magnetic pp col- 
lider detectors [0,0 , and the small tracking volume reduces backgrounds from muons from 



inflight decays of 7r and K mesons. 

The muon momentum p is measured from the muon deflection angle in the magnetic field 
of the toroid. The momentum resolution is limited by multiple scattering in the traversed 
material, knowledge of the magnetic field, and measurement of the deflection angle. The 
resolution in 1/p is approximately Gaussian and given by cr(l/p) = 0.18(p — 2)/p 2 © 0.008 
(with p in GeV) for the algorithm that was used to select the data presented here. The 
first of the two components in the above resolution function arises from multiple-Coulomb 
scattering in the iron toroids, and is the dominant effect for low-momentum muons. The 
second component is from the resolution on the measurement of the muon trajectory. 
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III. PARTICLE IDENTIFICATION AND EVENT SELECTION 



Because it is more difficult to separate the hadronic decays of W and Z bosons from 
the large background of dijet production, the cross section analysis uses the leptonic decay 
modes, W — > and Z — > 11 with I = e fl5fl , fi [[U|. 

The leptonic decays are characterized by a high-p^ lepton and large Ifix or by two high- 
ly leptons, for W or Z boson decays, respectively. This section describes the identification 
criteria used for electrons, muons, and neutrinos in this analysis. 



A. Events with Electrons 

Electrons are identified primarily by the presence of an electromagnetic shower in the 
calorimeter. A clustering algorithm finds these showers, and quality criteria are used to pick 
out electrons and photons, and thereby reduce backgrounds. Information from the central 
tracking system is used to separate electrons from photons. 

This analysis considers electrons in the Central Calorimeter (CC) defined by < 1.1, 
and the End Calorimeters (EC), 1.5 < \r)\ < 2.5. The region between the calorimeters 
is excluded because of poor resolution. In the CC, we also exclude electrons within 0.01 
radians in <p of the crack between adjacent calorimeter modules. 



1. Clustering Algorithm 

A nearest-neighbor cluster-finding algorithm [pjj is employed to find the electromagnetic 
energy clusters to be associated with electrons or photons. In each At] x A<ft = 0.1x0.1 tower, 
we sum the energies in all layers of the EM calorimeters. We then loop over all such EM 
energy towers with E > 50 MeV, and search the nearest-neighbor towers for high transverse 
energy. If there are other towers with E > 50 MeV, a local connection is made between those 
neighboring towers. In the next step, clusters are defined as groups of connected towers. If 
the transverse energy of the cluster is greater than 1.5 GeV, the cluster is saved for further 
analysis. The energy in the EM portion of the calorimeter is also required to exceed 90% of 
the total energy of the cluster, and the energy outside the central tower must be less than 
60% of the total. Both of these requirements are chosen to select clusters corresponding to 
narrow EM particle showers, as expected for electrons or photons. 

At this stage, the "electron" sample has a very large background from QCD processes 
(such as dijet production). This is because hadronic showers from QCD jets can sometimes 
fluctuate to look like electron or photon showers. Several other variables are introduced 
therefore to clean up the electron and photon selections. These variables involve a x 2 f° r the 
shape of the shower, the shower's isolation, and the spatial match between the calorimetric 
shower and the extrapolated position of some charged track emanating from the interaction 
vertex. 
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2. Covariance Matrix for the Shower 



Development of electron or photon showers in calorimeters is characteristically different 
from that of jets. The profile of the shower both in the longitudinal and transverse directions 
can therefore be used to discriminate between signal and QCD background. A covariance 
matrix is constructed to compare the shape of the experimentally observed shower with 
that expected from electrons or photons, taking into account the correlations among energy 
depositions in all the calorimeter cells in the cluster. 

For a sample of N electrons, a covariance matrix is defined as 

1 N 

M ij = ITt Y.( X ™~ < X i >)( X nj~ < Xj >), (3.1) 

^ V 71=1 

where x n i is the value of the zth observable for the nth electron, and < x\ > is the sample 
mean for that observable. There are 41 variables in the matrix: the fractions of shower 
energy in EM layers 1, 2 and 4, the fraction of shower in each cell in a 6 x 6 array centered 
at the hottest tower in EM layer 3, the logarithm of the shower energy, and the z position 
of the event vertex. A matrix , based on Monte Carlo simulation of electron showers, is 
constructed for each of the 37 detector towers at different values of \rj\. The Monte Carlo 
simulation was tuned to agree with the shower shapes of test beam electrons. 
For showers in the data, we calculate the H-matrix x 2 function: 

X 2 = XX X ~ < X 'i >) H ij( x J- < x 3 >)> ( 3 - 2 ) 

y 

where Xi is the measured value of the ith observable, and H = M~ l . Figure |2] shows 
the distribution of x 2 f° r showers from electron candidates from Z — > ee decays, and EM 
clusters in inclusive jet events that are primarily from overlaps between charged and neutral 
particle and 7r° decays. The two distributions are clearly different. Note that the covariance 
X 2 parameter will not necessarily follow a standard \ 2 distribution, because, in general, 



the observables defining the matrix are not normally distributed ||18|| . We require that an 
acceptable electron shower have x 2 — 100. 



3. Isolation Parameter 

An isolation variable is very useful for discriminating between background from jets and 
electrons from W or Z decay. Such electrons usually have very few other particles in their 
vicinity, while a jet contains many collimated particles close to each other. We therefore 
reject electron candidates with a significant amount of energy deposition nearby in the 
calorimeter. 

The isolation parameter for a cluster is defined by the fraction of energy in the vicinity 
of the core towers of that cluster 

, £ tot (^ = 0.4)-E EM (7e = 0.2) 

/is ° " E EM (K = 0.2) [6 - 6) 

where E tot is the total energy in the calorimeter in a cone with a radius 1Z = 
J (Arj) 2 + (A0) 2 = 0.4, and E EM is the energy in the EM section in a cone with radius 
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FIG. 2. x 2 distribution for (a) electrons and (b) jets. The arrows indicate the position of the 
cutoff used in this analysis at \ 2 = 100. 



0.2. Figure |3| shows the distributions found for electrons and jets. We require that accept- 
able electrons satisfy /i so < 0.10. 



4- Track Matching in the Central Detector 

Track information is used to distinguish electrons from photons. A reconstructed track 
is required to be within a 0.1 x 0.1 cone pointing towards the centroid of the EM shower. 
If this requirement is satisfied, the cluster is classified as an electron candidate, otherwise 
it is considered as a photon candidate. A track significance is defined as a measure of the 
quality of the match between the track and the centroid of the shower. For the central drift 
chambers (CDC), it is defined as: 



0"trk 



\ 




(3.4) 



while for the forward drift chambers (FDC), it is defined as: 



°~trk 



(3.5) 



9 



8 35 



o 
o 

c 

> 



30 
25 
20 
15 
10 
5 








(a) Electrons 




n n I n n 



_□ i □ nJ 



0.1 



0.2 0.3 
Isolation f. 



0.4 



0.5 




0.2 0.3 
Isolation f. 



ISO 



0.5 



FIG. 3. Distribution in the isolation fraction for (a) electrons and (b) jets. The arrows indicate 
the cutoff used in this analysis at /; so = 0.1. 



where z, p, and <fi are cylindrical coordinates and all the differences (J's) are calculated 
between the extrapolated coordinates of the track and the centroid of the shower in EM layer 
3 of the calorimeter. The standard deviations (ex's) in the denominators are the experimental 
resolutions for the corresponding matching parameters. Figure f| shows the distributions in 
<r tr k for electrons and jets in the CC, and the cutoffs used to define electrons. 



5. Energy Scale Calibration 



The absolute EM energy scale of the D0 calorimeter is determined using events in 
the Z — > ee mass peak An initial calibration was performed based on test beam 

studies of electrons and pions in a prototype calorimeter module. These determined that 
any nonlinearity and energy offsets of the calorimeter were negligible. We set the absolute 
scale by measuring the Z invariant mass peak and scaling our initial result to the known 
value of the Z mass [^(J. This correction is determined separately for each of the three 
calorimeter cryostats. The magnitude of the correction ranges from 1% to 7% [21 
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FIG. 4. Significance track matching for (a) electrons in the CC and (b) jets in the CC. The 
arrows indicate the value of the cutoff used in this analysis at <r tr k = 5 in the CC. 



6. Defining Electron Categories 



We define three categories of reconstructed electrons, referred to as "tight" , "standard" , 
and "loose". The tight criteria are used to reduce backgrounds as much as possible, while 
the loose criteria are used to obtain a higher reconstruction efficiency for electrons. The 
standard criteria are employed in the measurement of electron efficiencies. 

Tight electrons are defined as reconstructed EM clusters that: 

• pass the single-electron trigger (see next section) 

• have large EM fractions: / EM = f=u > 0.95 

• have H-matrix \ 2 < 100 

• are isolated: / iso < 0.10 

• have a good matching track, with <7 tr k < 5 for a CC, and o" tr k < 10 for EC. 

Standard electrons are defined with the same criteria, except for relaxed requirements on 
electromagnetic fraction (/ EM < 0.9) and isolation (/ iso < 0.15). The loose electron definition 
is the same as that for tight electrons, with the omission of the requirements for the trigger 
and for a matched track. 
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7. Single- Electron Trigger 



D0 uses a multiple-level trigger system. Common to many triggers used in this analysis 
is the Level-0 trigger, which requires signals in two hodoscopes of scintillation counters that 
are mounted close to the beam region on the front surfaces of the end calorimeters. Each 
analysis uses its own subset of Level- 1 (hardware) and Level-2 (software) triggers. 

A single electron trigger is used for both W — > ev and Z — > ee events to benefit from 
cancelations in trigger efficiencies when the cross section ratio is determined. The Level- 1 
trigger for single electrons demands that there be at least one electromagnetic trigger tower 
with transverse energy above 10 GeV (or 12 GeV for a small fraction of the early data). A 
trigger tower consists of four fixed calorimeter towers, covering Arj x A<p = 0.2 x 0.2, and 
contains most of the energy of an EM shower. 

The Level-2 trigger for electrons searches for the tower (Arjx A<p = 0.1 x 0.1) that contains 
the highest energy in the calorimeter, and then uses the nine (3 x 3) towers centered on it 
to form a cluster. The transverse energy for this cluster is required to be greater than 20 
GeV in order to pass Level-2. 

Level-2 also has minimal quality cuts on the shower shape of the cluster. The fraction 
of the cluster energy in the EM section must be above a given threshold which is dependent 
on the energy of and the position of the cluster in the detector. The transverse shape 
classification is based on the energy deposition pattern in the third EM layer. The difference 
of the energy depositions in two regions, covering Arj x A0 = 0.25 x 0.25 and 0.15 x 0.15 
and centered on the cell with the highest Et, must be in a window, which depends on the 
total cluster energy. Additionally, there is an isolation requirement, similar to that described 
above. The size of the outer cone in the trigger was set to either 0.4 or 0.6, with roughly 
half the data taken under each condition. 

8. Criteria on Shower Quality and Electron Kinematics 

Both the W and Z selections require one tight electron as defined in the previous sec- 
tion. The W — > ev selection requires a tight electron with E T > 25 GeV, 1$t > 25 GeV, 
and no second high Et electron electron. A total of 10338 candidate events satisfy these 
requirements. 

For Z events, there is an additional requirement of a second loose electron with Et > 
25 GeV. Also, the invariant mass of the two electrons (M ee ) is restricted to the range 75-105 
GeV. A total of 775 Z candidates satisfy all the criteria. 

Distributions of the transverse mass for the W — > ev events and the invariant mass for 
the Z — > ee events are shown in Fig. [| The transverse mass, M T , is defined by Mf = 
2Ej,^!t(^ — cos Afyv), where A<pi v is the azimuthal separation between the charged lepton 
and the missing transverse energy vector. 

B. Events with Muons 

Muons are identified by reconstructing a track from hits in the muon PDTs. The track 
is confirmed using information from the calorimeter and the central detector. Part of the 
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confirmation is a "global fit", which uses not only hit positions from the PDTs but also 
those from the other detector systems. It consists of a fit to the position of the primary 
vertex, a track in the central detector, and the muon track before and after the toroid. 
The seven parameters in the fit include four for the position and angle of the track before 
the calorimeter (in both the bend and non-bend views), two describing the effects of the 
multiple scattering in the calorimeter; and one for the reciprocal of the momentum (1/p). 
The momentum is therefore determined by the deflection in the magnetized iron with a 



correction for the expected energy loss in the calorimeter [22 



This analysis uses muons contained entirely in the central WAMUS detector (\rj\ < 1.0). 
To obtain a reliable momentum measurement, the minimal value of the integral of the 
magnetic field along the muon track is required to be / Bdl > 2 Tm. Although this re- 
duces significantly the acceptance for muons, it also eliminates a potential background from 
punchthrough. In the regions of low / Bdl, the D0 detector has only about 9 interaction 
lengths, while it has typically 13-18 interaction lengths elsewhere. This requirement there- 
fore provides a good momentum measurement, and a cleaner sample of muons because of 
the greatly reduced probability of hadron punchthrough for tracks from the calorimeter. 



1. Confirmation from Calorimeter and Central Detector 

Candidate muon tracks found in the PDTs must be confirmed by the presence of energy 
deposited along their trajectories in the calorimeter. This reduces background from cosmic 
ray muons and from random combinations of PDT hits. We require a sum of at least 1.0 
GeV of energy deposition in the cells of the extrapolated trajectory of the muon and in the 
two nearest-neighbor cells. A muon typically deposits ~ 3 GeV in this volume. Figure |5| 
shows the energy deposited in the calorimeter for good muons and for background. 

Another effective way to reduce background is to require a track match between that 
of the muon system track and the central detector. We require that there be a CD track 
associated with the muon, and that the angles between the two tracks match to within 
(muon track, CD track) < 0.04 radians, and A8 (muon track, CD track) < 0.12 radians. 



2. Rejection of Cosmic Rays 

Additional rejection of cosmic-ray muons and background from combinations of random 
PDT hits is provided by requiring small impact parameters of the muon track relative to 
the interaction vertex, and correct drift time relative to beam crossing. 

The impact parameters for muon tracks, both in the bend and nonbend views, are 
calculated by extrapolating the muon trajectory inside the toroids back to the primary 
interaction vertex. To be acceptable, the muon track has to point to the primary interaction 
vertex within 15 cm in the bend view and 20 cm in the nonbend view. 

The muon timing is determined by allowing the drift times of all the muon hits to vary 
coherently. The time interval t Q is defined as the offset between the beam crossing time 
and the time that gives the best fit for the muon track. Because they are produced in 
coincidence with beam crossings, prompt muons have a t Q distribution that peaks at zero 
However, cosmic rays arrive at random times. To have most of its PDT hits recorded, a 
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FIG. 5. Calorimeter energy for (a) good muons and (b) background. The arrows indicate the 
cutoff used in this analysis at E ca \ = 1 GeV. 



muon has to arrive within about ±400 ns of the beam crossing time (the total PDT drift 
time is ~ 750 ns). Due to the finite rise time of the trigger signals, the probability for 
accepting cosmic rays is enhanced somewhat for early arrivals (t^ > 0). To reject cosmic 
rays we require t < 100 ns. 

Figure ^ shows the distribution of t Q for signal and background samples. The background 
was obtained by selecting collider events containing two isolated high pt muons that are 
back-to-back (A#(/ii,/! 2 ) > 170° and A0(yUi,yU 2 ) > 160°). Such a sample is dominated by 
cosmic-ray muons. 



3. Global Fit to Muons 

The quality of the global fit of a muon track is characterized by the value of the y 2 for the 
fit, and depends on the parameters of the muon system as well as on those of the tracking 
system. By using the additional information, we are able to reduce the backgrounds from 
cosmic rays and from random combinations of PDT hits. The x 2 distributions for a signal 
and background are plotted in Fig. |7|. To accept an event, we require a fit with \ 2 < 100. 
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FIG. 6. The tg distributions for (a) good muons and (b) background. The arrows indicate the 
cut used in this analysis at t Q = 100 ns. 



4- Muon- Isolation Parameters 



A background that is not affected by the above criteria is that from QCD jet production. 
These events can have muons resulting from semileptonic decay of produced hadrons (e.g., 
bb events). Such muons are usually associated with jets, while muons from W or Z decay 
are most often isolated. We reduce the QCD background by imposing specific requirements 
on the calorimeter energy deposited within 1Z = 0.2 and 0.6 of the muon. 

We define the variable X 2 as the difference between the calorimeter energy observed in 
cells traversed by the muon (including the two nearest-neighbor cells within 1Z = 0.2 of the 
muon) and the expected contribution from the muon ionization, divided by the uncertainty 
in the expected energy loss in the calorimeter: 

_ £ tot (ft = 0.2)-£(expected) , Q 
2 2 = • (3.6) 

C"i?(expected) 



The expected energy loss is determined from the GEANT [£2| simulation of the D0 detector. 
We also define the variable Xq as 

Z 6 = E tot (n = 0.6) - E tot (n = 0.2). (3.7) 
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FIG. 7. Global-fit x 2 distribution for (a) good muons and (b) background. The arrows indicate 
the cutoff used in this analysis at x 2 = 100. 

Figure § shows the distributions of X 2 and X 6 for samples of isolated and nonisolated 
muons. The isolated muons are from a subset of W — > \w candidates with no jets opposite 
the muon in and the nonisolated muons are from events with muons in the range 10 < 
Pt < 15 GeV, a sample dominated by heavy quark decay. We reduce the QCD background 
significantly by requiring that X 2 < 3 and Iq < 6 GeV. 

5. Definitions of Muon Quality 

We define a "tight" muon as a reconstructed track in the PDTs that has: 

1. calorimeter confirmation with energy in central and nearest-neighbor cells of > 1 GeV 

2. a track match in the central detector 

3. a successful global fit, with x 2 < 100 

4. isolation requirements X 2 < 3 and X 6 < 6 GeV 

5. no back-to-back muon tracks (or PDT hits). 

A "loose" muon is not required to satisfy criteria 2-5. 
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6. Single-Muon Trigger 

The single muon trigger requires a high-p^ WAMUS muon candidate at both Level 1 and 
Level 2. The muon Level 1 system has two sublevels of hardware. The first sublevel passes 
events if there are PDT hits within a wide road (~ 60 cm), equivalent to a cutoff of 5 
GeV. The second sublevel searches in narrower roads (~ 30 cm), equivalent to a cutoff 
of 7 GeV. The Level 2 software trigger has pattern recognition, and accepts muons passing 
Pt > 15 GeV. Loose quality criteria are also applied at Level 2. 

Cosmic ray muons are suppressed at Level 2 if there is evidence of a single muon pene- 
trating the entire detector. Muon candidates with a track in the opposite muon chambers 
within 20° in and 10° in 9 are rejected, as are those candidates with PDT chamber hits 
on the opposite side within 60 cm (roughly 5°) of the projected muon track. 

7. Muon Kinematic and Quality Criteria 

The W — > p,v offline selection requires one tight muon with pr > 20 GeV and Ifr > 
20 GeV. The Z — > fifi offline selection requires at least one tight and one loose muon. Both 
muons must have px > 15 GeV and at least one has to have pt > 20 GeV. To reject the 
cosmic ray background, we require either A0 < 160° or A8 < 170° between the two muons. 
To eliminate low-mass dimuon pairs it is required that A<j) > 30°. Any event satisfying the 
Z criteria is removed from the W sample. 

Distributions of the transverse mass for W events and the dimuon invariant mass for Z 
events are shown in Fig. ||. 

C. Neutrino Identification 

Neutrinos are identified in the D0 detector by the presence 
(#r). We define 

where 



^T y = -^2 Ei sin 6i sin 4>i (3.10) 

i 

where i runs over all calorimeter cells with readout signals after zero suppression, and Ei 
is the energy deposited in the ith cell, with d{ and as the polar and azimuthal angles of 
that cell, respectively. If there are muons in the event, we subtract the px of the muons as 
follows: 

$T X = - E i sin 0t COS <f>i - p Tx , (3.11) 



of missing transverse energy 

(3.8) 
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FIG. 9. Transverse mass distributions in the final (a) W — > ev and (b) W fj,v samples and 
invariant mass distributions in the final (c) Z — > ee and (d) Z — ► ^u/x samples. The points are the 
data, the hatched regions correspond to the estimated backgrounds, and the histograms are the 
sum of signal (Monte Carlo) and background. 
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$ T y = ~J2 E i sin 9i sin <t>i ~ PTy (3- 12) 

i 

The resolution of the missing transverse energy is affected by many factors, such as 
statistical energy fluctuation in the calorimeters, energy lost in and around the beam pipe 
and cracks in the central calorimeter, signal fluctuations caused by the uranium radioactivity, 
random and coherent electronic noise. 

Since we have a nearly hermetic calorimeter with good energy resolution, we also obtain 
very good J^t resolution. A global quantity called the scalar transverse energy, defined as 

Y / E T = Y, E i sin9 i ( 3 - 13 ) 

i 

is used to parameterize the resolution as 

^ T =a + &5> T (3.14) 



with a = 1.08 GeV and b = 0.019 obtained from minimum-bias data. Figure |10| shows the 
dependence of the 1£t resolution on Et- 




IV. BACKGROUNDS 

Backgrounds to W and Z events can be into two groups: those from "fake" leptons, 
whose levels are estimated from data, and those from "physics" processes that contain true 
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isolated high-p^ leptons and true I$t- The contributions for the latter sources are estimated 
from Monte Carlo samples. 

Electron background stems primarily from jets and direct photons passing our electron 
criteria. Muon background consists mainly of cosmic-ray muons, random hits in the muon 
chambers that form a track, and muons from heavy quark decays. The inherent background 
processes, common to both lepton channels, are W — > rv — ► Ivvv, Z — > 11, Z — > tt —>■ llvv, 
and Drell-Yan production of pairs. 

A. Backgrounds to W — > Iv 

1. QCD background to W ^ ev 

A multijet event can be misinterpreted as a W — > ev in our detector, for instance, if one 
of the jets fluctuates to have a high electromagnetic content and passes the electron selection 
requirements, while another jet loses energy in the cracks of the detector, or its energy is 
otherwise mismeasured, to yield missing transverse energy, thereby faking a neutrino. 

We study this background through the $r distribution of tight-quality electrons prior 
to the imposition of any If? criteria. Figure [11] shows that there are two peaks in the data. 
The peak in the low-$r region is mostly due to jet events, and the second peak is dominated 
by true W — > ev decays. 




10 20 30 40 50 60 

Missing Transverse Energy (GeV) 



FIG. 11. The rfx distribution for the sample of W — > ev events (solid histogram). The dotted 
points with error bars represent the multijet QCD background normalized in the low If? region 
(0-10 GeV) to the signal sample. The arrow indicates the cutoff for the final selection of I$t > 25 
GeV. 
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We also consider a sample of QCD dijet events, those for which the electron candidate 
fails the isolation criterion (i.e., we require that there be some energy deposition around the 
"electron", which is presumably due to the rest of the remnants of the jet). Since isolation 
and I$t criteria are not correlated, the I^t spectrum for the dijet events in this sample and 
in the tight-electron sample should be the same. We therefore normalize the two samples 
in the low-^r region, and extrapolate to find the number of background events under the 
W peak passing the IpT cutoff of 25 GeV. Figure [II] shows the second background sample 
normalized to the tight-electron sample. Since the background falls rapidly with , there 
are very few events that pass the cutoff. 

We consider the QCD background to W decays separately for electrons found in the 
CC and EC calorimeter. The events are further subdivided into two groups, to take into 
account two variants used in the electron trigger (with isolation radius 0.4 and 0.6). For 
each of these data subsets, the background sample is normalized to the signal sample in 
the region 0-10 GeV in If?, to determine the number of events in the background samples 
with If? > 25 GeV. We find the background to electrons in the CC to be ~ 3% and in the 
EC to be ~ 4%. The overall jet background in the W — » ev data sample is (3.3 ± 1.7)%. 
The uncertainty has changed with respect to the original Letter || as a result of additional 
studies of the $r distribution pi 

2. QCD background to W — > ixv 

The QCD background to W — > fiu events consists of muons from decays of particles 
associated with jets. Most such muons fail our isolation criteria. We estimate the background 
by fitting the observed distribution in energy deposition X 6 (without imposing any isolation 
criteria) to a sum of distributions expected for isolated and for nonisolated muons (see 
Section III.B.4). The fit to a linear sum of signal and background to the data is shown 
in Fig. [12|. After applying the two isolation criteria, the QCD contamination in the final 
W ^ \iv sample is (5.1 ± 0.8)%. 

3. Backgrounds toW^^LU from Cosmic-Ray Muons and Random PDT Hits 

The background from cosmic rays and random PDT hits is estimated from the t Q dis- 
tributions. Since neither background is beam-associated, there should be no correlation 
between the best time for the fit of the track and the beam crossing. The prompt distribu- 
tion is obtained from a sample of muons with > 5 GeV, and very tight quality criteria: 
a matching track in the central detector, tight global fit x 2 , and sufficient energy deposition 
around the muon trajectory to ensure that it is part of a jet. The background fraction is 
determined by fitting the data sample to a linear sum of signal and background v Q distribu- 
tions. The contamination from cosmic rays and random hits in the final W — > \xv sample is 
estimated to be (3.8 ± 1.6)%. 
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FIG. 12. The distribution in the isolation X% for the W — ► fiv sample (data points), fitted to a 
linear sum of the signal and background. The dashed histogram corresponds to isolated muons, the 
dotted histogram shows muons originating from QCD processes, and the solid histogram represents 
their sum. 

4- Punchthrough and tt/K Backgrounds to W —* fiu 

For pt > 12 GeV, the background originating from ir/K decays is estimated to be an 
order of magnitude smaller than that from b decays |[24|| . The rate from punchthrough is 
expected to be yet another order of magnitude lower. The reason for the low rate is the 
great thickness of the calorimeter and iron toroid systems at D0, combined with the fact 
that the momentum measurement is made after most of the material has been traversed. 
The background contamination from these sources in the W — > \w and Z — > samples is 
therefore negligible. 

5. W — > tv — ► Ivvv Backgrounds toW^lv 

The process W — > tv — > Ivvv is experimentally indistinguishable from the signal. There- 
fore the only means for reducing this background is through differences in kinematics. Since 
the background (charged) lepton comes from the decay of a r, it will have a much softer 
distribution than from direct W decay. Just the standard kinematic requirements keep this 
background to a moderate level. 

We use Monte Carlo simulations to calculate the geometric and kinematic acceptance A T 
of W — > tv — > e([i)vvv. Accounting for the r — > e(/x) branching fraction, we find the overall 
r background in the electron channel to be (1.8 ± 0.2)%, and in the muon channel, it is 
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(5.9 ±0.5)%. The background in the electron channel is lower because the energy resolution 
for electrons is better than for muons and the px cutoffs are higher in the electron analysis. 

6. Z — ► 11 Backgrounds toW^lv 

One of the two leptons from Z decay can escape detection or be poorly reconstructed in 
the detector and thereby simulate the presence of a neutrino, and contribute to the W — > Iv 
data sample. Assuming a{pp — > Z — > ee)/a(pp — > VF — > ez/) is 0.10 ± 0.01, we find from 
a Monte Carlo simulation that this background fraction is (0.6 ± 0.1)% for electrons, and 
(6.5 ±0.5)% for muons. The electron background is lower because of the greater hermeticity 
of the calorimeter for electrons compared to WAMUS for muons. 

7. Z — ► tt — > Ivvlvv Backgrounds toW^lv 

The process Z — > rr has the same rate as Z — > //, which is already ten times smaller than 
the rate of production. Each electron from r decay has the soft pt spectrum mentioned 
for the case of W — > tza This background is therefore doubly suppressed. For the muon 
channel we estimate the background to be (0.8 ± 0.2)%, and for the electron channel it is 
negligible. 

B. Backgrounds in the Z — > 11 Sample 

1. QCD background to Z — > ee 

The background to Z — > ee consists mainly of QCD jet production, where the jets are 
misidentified as electrons. Because the invariant mass distribution of two electrons from the 
Z decay has a well defined resonance peak, and the background has only a weak dependence 
on the ee mass, we use the shape of the mass distributions to estimate background. We fit a 
theoretical Z/7 line shape and the experimentally determined shape of the QCD background 
(see below) to the data, and determine the absolute normalization of the QCD background 
through this fit. 

The invariant mass distribution of the QCD background is obtained from data. We can 
approximate the "two-electron" mass spectrum from jets with p? > 25 GeV, for the mass 
range 65-250 GeV, by an exponential function: 

/ dijet (m) oc e -°- 0237m . (4.1) 

A second QCD contribution arises from direct-photon events with associated jets. Here 
the jet fragmentation fluctuates sufficiently for the jet to be reconstructed as an electron, 
while the photon is mistaken as a loose electron (only failing the track match). Again we 
can describe this "dilepton" mass spectrum by the following exponential function for the 
same mass range as above: 

/ 7jet (m) oc e ^.034 5m _ (4 . 2 ) 
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We use the PYTHIA Monte Carlo |2S ] to generate the complete Zj 7 line shape, including 
QED radiation from electrons, and we also simulate the energy resolution of the detector. 

Using a maximum likelihood fit to the dielectron invariant mass spectrum, we determine 
the fraction of events in the data sample that can be attributed to QCD background. We find 
a total background of (2.8 ± 1.4)%, where the error includes statistical as well as systematic 
uncertainty to account for the sensitivity to the mass window used in the fit (71 or 75 GeV 
to 111 or 121 GeV). 



2. QCD Background to Z — > /i/i 



For the Z — > /i/j, sample, the background is estimated in a similar fashion to that for 
the W — > iw background, by fitting the calorimeter energy distribution X 6 . The QCD 
background in the final Z — > \iyu sample is estimated to be (2.6 ± 0.8)%. 



3. Backgrounds to Z — > /j/i from Cosmic-Ray Muons and Random PDT Hits 



The backgrounds from cosmic rays or random hits in Z — ► /i/x are estimated from the 
muon track-time distributions (tl) using the same fitting techniques as used for the W- 
background estimate. The total contamination from these sources in the final sample is 
found to be (5.1 ±3.6)%. 



4- Z — ► tt Background to Z — > 11 



The process Z — > tt, where both taus decay to either electrons or muons is a small 
background in this analysis. In the Z — > ee sample, the reduced acceptance (due to the 
soft pt spectra of the electrons) and small r branching fractions allow us to neglect this 
background. For the Z — > nn sample, we estimate a background of (0.7 ± 0.2)%. 



5. Drell- Yan Pair Production 



The Drell- Yan process 7 — > 11 is coherent with Z — > 11, and the experimentally observed 
production of lepton pairs corresponds to the square of the sum of the M 7 and Mz ampli- 
tudes. However, because we are interested in comparing to theoretical predictions for the 
\Mz\ 2 term, the size of the Drell- Yan fraction |M 7 | 2 and the interference term must therefore 
be deduced before making our comparison. 

We use the ISAJET Monte Carlo [26] to estimate these two terms relative to pure Z 
production, and express them as the fraction of the number of pure Z events. We cross 
check with PYTHIA [p5l and find a similar number. This "background" is thus found to be 



[1.2 ± 0.1)% for the electron channel, and (1.7 ± 0.3)% for the muon channel. 
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Backgrounds 


W -> ei/ 


W ^ jJLU 


Z — ► ee 


Z — > /i/i 


QCD Dijets 


3.3 ±1.7 


5.1 ±0.8 


2.8 ± 1.4 


2.6 ±0.8 


Cosmic Kays, etc. 




3.8 ± 1.6 




5.1 ± 3.6 


W -> r -» e(/i) 


1.8 ±0.2 


5.9 ±0.5 






Z — > ee(fifi) 


0.6 ±0.1 


6.5 ±0.5 






Z — > TT — ► 




0.8 ±0.2 




0.7 ±0.2 


Drell-Yan 






1.2 ±0.1 


1.7 ±0.3 


Total 


5.7 ± 1.7 


22.1 ± 1.9 


4.0 ± 1.4 


10.1 ±3.7 



TABLE I. Backgrounds to the W and Z samples. All estimates are in per cent of the total 
candidate samples. 



C. Summary of Backgrounds in the W — > Iv and Z —> 11 Samples. 

The backgrounds are summarized in Table |. We find a total background in the W — > Iv 
samples of (5.7 ± 1.7)% for electrons and (22.1 ± 1.9)% for muons. For Z — > 11, the total 
background estimates are (4.0±1.4)% and (10.1 ±3.7)%, for the electron and muon channels, 
respectively. 



V. DETECTOR SIMULATION AND ACCEPTANCE 

The acceptances for the processes pp — > W — > Iv and pp — > Z — > 11 are defined as the 
fractions of all the W —>■ Iv or Z — > 11 events that pass our fiducial and kinematic criteria. 
The acceptance is estimated using an event generator to model vector-boson production and 
decay, and a Monte Carlo simulation of the D0 detector. This section describes the event 
generators, the detector simulations, and the results of the acceptance calculations. 



A. W — > ev and Z — > ee Simulation 

A fast Monte Carlo simulation is used to calculate the acceptance in the electron channel. 
A parton-level generator produces a vector boson, which is made to decay to leptons in the 
boson rest frame. The leptons are boosted to the lab frame according to the longitudinal and 
transverse momentum of the boson. The longitudinal momentum of the boson is determined 
by the parton distribution functions and the energy in the center of mass. The transverse 
motion is caused by radiation of initial-state partons or through higher order contributions 
toWorZ production. We calculate the pr spectra from the double differential cross section 



a ra/dprdij provided in a next-to-leading order (NLO) program |27j]. The calculation uses 
a standard perturbative method for high p?, a resummation scheme for the low pt region, 
and a matching scheme between the two. 

The specific procedure involves, first, generation of the rapidity of the vector boson from 
the randomly-selected momenta of the incident quarks. Then, the double differential cross 
section at that rapidity value is used to generate a px distribution, from which the px of the 
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vector boson is chosen. Once all the four-vectors of the W or Z boson and decay leptons 
are generated, the differential cross section is calculated and used as a weight for this event. 
One million such events were generated, and the weights used to obtain the acceptance of 
our geometrical and kinematical criteria. 

The detector simulation includes modeling of the primary vertex distributions, the elec- 
tron energy and Ifir resolutions, and the turn-on of the Level 2 trigger. The vertex (im- 
position) of collisions is generated from a Gaussian distribution with o = 30 cm and 
< z >= —8 cm, to reproduce the measured vertex distribution in the data. The electron 
energy is smeared with a resolution of 

+ (A)' 

where S = 15.7%VGeV and C is 0.4%. 

The missing E? is reconstructed from Ifi? = Pr ~ Pti where p^ is smeared according to 
the EM energy resolution, and the pt of the W is smeared to match the hadronic energy 
resolution, because p^ is determined from the hadrons recoiling against a W boson. A 
correction factor of 0.83 is applied to the hadronic energy scale of the calorimeter; this 
factor is obtained by studying the balance of the sum of the px of two electrons and the 
recoil hadrons in Z events. 

The effect of underlying events in the data is included in the detector simulation. A 
vector of chosen at random from a sample of minimum bias events, is added to the 
above Ifx to simulate the smearing contributed by the underlying event. 

The event simulation includes radiative corrections. We first calculate the acceptance for 
the radiative final states W'y or Z^ (with a threshold energy for the photon of greater than 
20 MeV) p8| . In this calculation, we define a cone whose axis is centered along the direction 
of the electron, with the cone size defined as 1Z = y/5r]' 2 + S(p 2 , where 5i] and 5(p are the 
differences in rj and <fi between the electron and photon directions. If R < 0.3, we add the 
energies of the electron and photon, and treat the sum as the electron energy. Otherwise the 
energy of the electron is left intact. We then combine the acceptances of W/Z and W'y/Z^f 
to get the final acceptances. The radiative corrections are 0.6% for our W measurement and 
1.6% for the Z channel. 

The final effect included in the simulation is the electron px cutoff at Level 2. The effect 
of this cutoff is folded in using the trigger turn-on curves observed in the data (see Fig. |T3"| ). 



B. W — * /iu and Z — > fifi Simulation 

A full detector simulation is used for the acceptance calculations in the muon channel. 
We use ISAJET [^(J as the parton-level event generator, and GEANT |22| to model the D0 



detector. We incorporate the PDT efficiencies and resolutions on a chamber-by-chamber 
basis using measurements from the data. 

The chamber efficiencies are obtained from large samples of good quality muons, and 
include variations observed during the run. We check these efficiencies by comparing the 
Level 1 trigger efficiencies predicted from the simulation with those observed in the data. 
The individual chamber resolutions are obtained from residuals in fits to muon tracks in 
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FIG. 13. Trigger efficiency for electrons in the CC as a function of reconstructed Ex- The 
threshold behavior is similar for electrons in the EC. The line is a function fitted to the data 
points. 



collider data. Finally, the overall momentum resolution in the Monte Carlo is tuned to fit 
the shape of the reconstructed Z — > and W — > fiis mass distributions. 



C. Acceptance Calculation 

We define the acceptance as the fraction of events passing our kinematic and fiducial 
requirements. For Z — * ee, we also include the effect of requiring the invariant mass of the 
lepton pairs to be within the mass window from 75 GeV to 105 GeV. This mass requirement 
accepts (95.5±0.3)% of the events. For Z — > fift, we also include the effect of the requirement 
on the angles between the two muons. This has an efficiency of only (81 ± 1)%, but greatly 
reduces the cosmic ray background. 

The overall acceptances A w for W — > lis, and A z for Z — > 11, are summarized in Table [TI|. 
The acceptance in the muon channel is much lower than for the electron channel because, 
for this analysis, we restrict the r] range to the region where the single muon trigger is most 
efficient. The effect is roughly a factor of two per lepton. 



D. Systematic Errors on the Acceptances 

We summarize the systematic uncertainties for acceptances of W and Z bosons in the 



electron channel, and for the ratio of acceptances for the two processes, in Table QL The 
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A w 


A z 


Electron Channel 


(46.0 ± 0.6)% 


(36.3 ±0.4)% 


Muon Channel 


(24.8 ±0.7)% 


(6.5 ±0.5)% 



TABLE II. W and Z acceptances in the electron and muon decay channels. The systematic 
errors come from a variety of sources, as explained in subsection V.D. 

errors on the ratio are calculated separately, taking account of the partial cancelation of the 
systematic errors in A w and A z . 

The systematics from the parton distribution functions (pdf) are studied by compar- 
ing results obtained with CTEQ2M, CTEQ2MS, GRV, MRSD0' and MRSD-. We define 
CTEQ2M as the central value of our calculations, recalculate the acceptances for each pdf, 
and quote the maximum difference in our prediction as the systematic error. We also vary 
the parameterization of the pt spectrum used to generate the W and Z events within the 
range consistent with our data. 

The values of the W mass and width |J are varied by one standard deviation and 
the consequences propagated through the acceptance calculations. The W acceptance is 
sensitive to the mass of the W boson, varying by 0.7% for a change of mass of 0.18 GeV. 
The result is not sensitive to the width of the W boson. The errors from uncertainties in 
the mass or width of the Z are extremely small, and we therefore neglect them. The error 
in the simulation of the Ifi? depends mainly on the smearing and simulation of the energies 
of soft jets. This error is dominated by the uncertainty in the hadronic energy scale. 

We vary the parameterizations of the Level-2 trigger efficiency and the input-vertex 
distribution to estimate the uncertainties from these sources. 

The electromagnetic energy scale error and resolution also contribute to the systematic 
uncertainty on the acceptances. We vary both of these within their measured uncertainties, 
and find that the uncertainties from both these sources are small. 

The error in the acceptances due to radiative corrections is obtained by varying the cone 
size used to estimate the corrections. 

The muon channels have quite different, and, in general, larger systematic uncertainties 



than the electron channels. These uncertainties are summarized in Table [TV. 

We determine the systematic error from the chamber efficiencies by varying the efficien- 
cies within their uncertainties and repeating the acceptance calculation. We determine the 
uncertainty in resolution by varying the overall resolution within the constraint that it be 
consistent with the observed Z — > \m mass distribution. 

The dependence on pdf is larger for the muon channel than for the electron channel 
because of the difference in r\ coverage. The muon rapidity range ends in a region of large 
acceptance, and so is quite sensitive to the input parton distributions, while the electron 
coverage is more complete, extending in pseudorapidity to a region where the contribution 



to the cross section is relatively small, and so is less sensitive to the inputs. Figure ITl shows 



the generated pseudorapidity distribution for W — > £u and accepted regions for muons and 
electrons. 
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A \Al -*f>1J 1 A 7, >Pf> 

/I V V ' L- IS I j\£j ' L 


Choice of pdf's 


0.4% 


0.6% 


0.3% 


W/Z pt spectra 


0.3% 


0.2% 


0.4% 


W mass 


0.7% 




0.7% 


W width 


< 0.2% 




< 0.2% 




0.6% 




0.6% 


Trigger efficiency 


0.3% 


< 0.1% 


0.3% 


Vertex distribution 


0.4% 


0.4% 


0.3% 


EM energy resolution 


0.1% 


0.3% 


0.2% 


EM energy scale 


0.3% 


0.3% 


0.2% 


Radiative corrections 


0.3% 


0.4% 


0.4% 


Z mass window 




0.3% 


0.3% 


Total 


1.3% 


1.0% 


1.3% 



TABLE III. Relative uncertainties in the acceptances for W 
of these estimates are given in the text. 



ev and Z — > ee events. Details 











Choice of pdf's 


2.0% 


4.0% 


2.0% 


Chamber efficiencies 


0.7% 


4.4% 


3.6% 


Chamber resolutions 


1.0% 


0.4% 


1.4% 


Monte Carlo statistics 


1.7% 


3.9% 


4.2% 


Total 


2.9% 


7.1% 


6.0% 



TABLE IV. Systematic uncertainties in the acceptance for W — > fxu and Z — > fjbfi events. All 



other systematic uncertainties are negligible compared to those listed in this table. 
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FIG. 14. Pseudorapidity distributions for charged leptons in W — > Iv events. The solid line is 
the generated distribution for either electrons or muons. The dotted and dashed lines show the 
distributions for the electron and muon channels, respectively, after applying the fiducial require- 
ments. 

VI. CALCULATIONS OF EFFICIENCY 

To measure the efficiency for particle identification (ID) and triggering requires choosing 
clean and unbiased samples of electrons and muons. We use the Z — > ee and Z — > fx/i 
samples as our source of high-p^ leptons with low background. Requiring only one lepton 
to pass all the particle ID and trigger criteria, leaves the other lepton unbiased with respect 
to these cutoffs, and it can be used to measure the efficiencies. 

A. Electron Efficiencies 

For obtaining electron efficiencies, we use the Z — >■ ee events in the peak region 86 < 
M ee < 96 GeV. As usual, we determine the systematic error by varying the electron selection 
criteria and the procedure used for background subtraction. We used two sets of selection 
cuts, one being the tight cutoffs, and one being the standard cutoffs (see Sect III. A. 6). 

The major background in the Z sample is from QCD jet production with two jets 
misidentified as electrons. The two-body invariant mass distributions are used to study 
this background. As discussed earlier, over a limited mass range, this background can be 
approximated by an exponential dependence, or by simple polynomials. We extract the 
background fraction in three ways: 
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• Method 1: The average of the number of events in the two sideband regions of the Z 
peak, 61 < M ee < 71 GeV and 111 < M ee < 121 GeV, is taken as the background in 
the peak region. 

• Method 2: The mass spectrum is fitted to the sum of a relativistic Breit-Wigner shape 
for the Z convoluted with a Gaussian resolution for M ee , and a linear function for the 
background. The result of the fit in the region 60 < M ee < 70 GeV is used to estimate 
the background. 

• Method 3: The same fit in Method 2 is used but the two sideband regions M ee < 70 
GeV and M ee > 110 GeV are used to measure the background under the Z peak. 

The efficiency of a particular set of criteria is measured by 

~ £ bfb la i \ 

£=—, j- I 6 - 1 ) 

1 - Jb 

where e s is the efficiency measured for events in the peak region, e b is that measured for 
events in the background region and /& is the fraction of background in the peak region. 
Note that it does not matter that there are some real Z — > ee events in the background 
region. Our measurement of the efficiencies is correct as long as there is only signal left in 
the peak region after background subtraction. 



1. Efficiency for Electron Selection 
The single-electron selection efficiency can be expressed as 

£ ^e_shwr ' ^e_trk ' ^e_trigj (6-2) 

where £ e ^hwr is the efficiency from the requirements on shower shape in offline reconstruction, 
£ e _trk is the efficiency of associating a good track with an electron cluster in the calorimeter, 
and £ e _tri g is the efficiency of the trigger. 

The £ e ^hwr values reflect the following selections 

• H-Matrix x 2 < 100 

• /iso < 0.10 

• Jem > 0.95 

Efficiencies for electrons are measured separately in the CC and in the EC. We extract six 
values, corresponding to the two sets of tagging criteria (standard and tight), and the three 
background estimations discussed above. The central value is taken to be the median value 
closest to the mean, and its uncertainty as the combination of the statistical uncertainty 
associated with the central value and half the difference between the highest and lowest 
calculated efficiencies. The results are 

£shwr_cc = 0.881 ± 0.015 



32 



and 



^shwr^c = 0.884 ± 0.024, 

where the errors are dominated by statistics. 

The efficiency for reconstructing a track associated with an electron has two components: 
the efficiency for finding a track near the electron, and the efficiency of the cutoff on the 
parameter <7 tr k, which provides the quality of the match between the position of the track 
the calorimeter shower. The efficiency of track matching is measured by taking the ratio of 
the number of calorimeter clusters in the Z — > ee sample that are reconstructed as electrons 
to the number reconstructed as either electrons or photons. (The only difference is the 
presence or absence of a track match.) 

The (j trk and trigger efficiencies are obtained in the same way as the efficiencies for 
shower shape cuts in the calorimeter. Combining e e _ tr k and e e _ trig , using our six estimations, 
we obtain 

£trk_cc • e trig _cc = 0.830 ± 0.014 

and 

£trk_EC • £trig_EC = 0.774 ± 0.024. 
The uncertainties are again dominated by the statistics of the Z event sample. 

2. Efficiencies for W — > ev and Z — > ee Events 

The selection criteria for electrons from W events are just the ones given above (tight). 
The total selection efficiency is therefore the product of the shower, track, and trigger terms, 
giving 

e^* ev = 0.731 ± 0.018 (6.3) 

and 

e^* ev = 0.684 ± 0.028. (6.4) 

We combine the CC and EC results, by weighting them by their relative acceptances, to 
obtain 

e w ^ eu = 0.704 ± 0.017. (6.5) 

For Z — > ee events, one of the two electrons is selected exactly as in W — > ev events. 
The second one is selected without imposing any track or trigger requirements, but with all 
the other criteria being the same; thus the total efficiency for selecting this loose electron is 

JUSt £shwr- 

The Z — > ee efficiencies are therefore 

e^c-cc = 0- 753 ± O-O 20 , 
eg£f| c = 0.748 ± 0.023, 
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and 



£ f^_ e | c = 0.742 ± 0.042. 

The overall efficiency is obtained by weighting using the relative acceptances of CC and EC 
events, giving 

£ z^ee = g 735 ± Q.024. (6.6) 

The ratio of efficiencies of W to Z selections is calculated directly using each of our 
different methods, and the systematic error assigned independently from e w or e z . In this 
way, any correlation of systematic errors of e w ^ ev and £ z ^ ee is taken into account. The 
ratio of efficiencies is 



£ Z-»ee 
£ W^eu 



1.045 ±0.019. (6.7) 



B. Muon Efficiencies 

The single-muon efficiency can be written as 

where e^_ reco is the muon reconstruction efficiency, e^m is the muon ID efficiency, and £ M _t r i g 
is the muon trigger efficiency. Each of these efficiencies is measured using a different unbiased 
data sample, as described below. 



1. Muon Reconstruction Efficiency 

The muon reconstruction efficiency is estimated using events from a special data run 
that had no Level-2 requirements. We require that there be a jet reconstructed offline in 
the same r]-<j) region as a muon candidate found by the Level- 1 trigger. No Level-2 or muon 
reconstruction criteria are imposed either online or offline, and the muon candidates are 
categorized as "good" or "bad" tracks through visual examination of the event displays. The 
efficiency is defined to be the percentage of "good" muon tracks which pass the loose muon 
ID (defined in Sec. |I 1 1 B 5|) . To reduce the systematic error of this method, the scanning is 



performed by at least two physicists. The offline reconstruction efficiency for "good" tracks 
is found to be 

£ M _reco = 0.952 ± 0.033 (6.9) 

per muon. 



34 



2. Muon Identification Efficiency 

Muon efficiencies are derived directly from the Z — > fifi sample. One muon tags the 
event by passing all the particle ID requirements, leaving the other muon to form part of 
an unbiased sample of isolated high-p T muons. The p T cutoff is raised to 20 GeV on both 
muons, which minimizes the backgrounds from QCD, cosmic rays, and random PDT hits. 

The overall muon ID efficiency is found to be 

^_id = 0.626 ± 0.047. (6.10) 

The most important contributions to the inefficiency stem from the requirement on track 
matching in the CD (e = 0.82 ± 0.04), for finding a track and matching the angles and from 
the calorimeter isolation criteria (e = 0.85 ± 0.04). 



3. Muon Trigger Efficiency 

The muon trigger efficiency is not estimated from the Z — > sample because of poor 
statistics. Instead, a sample of "unbiased" muons that pass quality criteria, and are present 
in events passing a non-muon (usually a jet) trigger, are used. The results are cross-checked 
with the ones obtained from the Z — > sample, and the results agree within the statistical 
uncertainty. 

The overall trigger efficiency for high-p T muons is 

^-trig = 0.367 ± 0.019. (6.11) 

The relatively limited geometric coverage of the muon chambers is the most important factor 
contributing to this low efficiency. The Level- 1 trigger requires hits in all three layers of the 
PDT system, and only m 60% of muon tracks in the fiducial region satisfy this requirement. 
For tracks that satisfy it, the trigger efficiencies are (80 ± 2)% for Level 1 and (78 ± 3)% for 
Level 2. 



4- W — > [xv and Z — ► fjbfi Efficiencies 

The efficiency for W — > [iv candidates corresponds to the single-muon efficiency described 
above. The total efficiency is therefore the product of the reconstruction, ID, and trigger 
terms, and equals 

£ w-+nv _ Q 219 ± 0.026. (6.12) 

The Z — > fifi efficiency takes into account the fact that both muons must be recon- 
structed, but that only one has to pass the ID and trigger criteria. Combining the single 
muon efficiencies, we obtain 

gZ-w* = 0.527 ±0.049. (6.13) 

The ratio of efficiencies ofWtoZ selections takes account of the correlations among the 
systematic errors. The result is 

— = 2.48 ±0.19. (6.14) 
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VII. CROSS SECTION 



A. Determination of Integrated Luminosity 



The integrated luminosity is determined by monitoring non-diffractive inelastic pp col- 
lisions using two hodoscopes of scintillation counters (the Level-0 trigger [T^|) mounted on 



the front surfaces of the end calorimeters near the beam axis. The average |29[ of the values 
measured by the CDF [30] and E710 |M| experiments at Fermilab is used for the inelastic 
cross section. The reaction rate measured by the Level-0 system corresponds to a cross 
section (the Level-0 visible cross section) of cr L0 = 46.7 mb. 

For the electron trigger used in this analysis, after corrections for experimental dead 
times and multiple interactions, the integrated luminosity is determined to be 



J Ldt = 12.8 ± 0.7 pb 
while the muon trigger had an exposure of 



Ldt = 11.4 ±0.6 pb" 



(7.1) 



(7.2) 



The 5.4% systematic uncertainty in the luminosities is calculated from the uncertainty in 
the pp inelastic cross section (4.6%), the systematic errors on the acceptance (2.0%), and 
efficiency (2.0%) of the Level-0 detectors. 



B. W and Z Production Cross Sections 

1. Theoretical Predictions 

The W and Z boson total production cross sections have been computed from a complete 
calculation to order a 2 s [0. We used M z = 91.19 GeV, M w = 80.23 GeV, and sin 2 9 W = 
1 - Myy/Mz = °- 226 ' and CTEQ2M @ pdf for our central value and considered the other 
pdf sets shown in Table |V[ The strong correlation of the W and Z boson cross sections 
decreases the sensitivity of the ratio of cross sections to variations of the pdf. Taking 
CTEQ2MS and CTEQ2ML as the extremes, we obtain a w /a z = 3.33 ± 0.02. 

Until recently, the uncertainties on the calculated cross sections were dominated com- 
pletely by the variation due to choice of pdf. Recent measurements of the proton structure 
function F 2 and of the W ± rapidity distributions have restricted the acceptable pdf choices 
to the point that other sources of error must be considered. The sources we considered are 
the use of NLO pdf sets instead of NNLO (which would be more appropriate for use with 
the 0{a 2 s ) calculation), variation of the calculated cross section from the uncertainty in Mw, 
and the uncertainty due to the dependence on renormalization and factorization scales. 

While the W and Z boson total cross sections have been calculated up to 0(a 2 ), the 
corresponding NNLO pdf sets are not yet available. The uncertainty on the cross sections 
due to using NLO pdf's has been estimated |34| to be 3% at y/s = 1.8 TeV. This uncertainty 



is assumed to cancel in the predicted ratio of cross sections. 
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pdf 


a w (nb) 


o z (nb) 




MRSSO/ 


22.114 


6.633 


3.334 


MRSDO/ 


22.150 


6.680 


3.316 


MRSD/- 


21.810 


6.558 


3.326 


MRSH 


22.043 


6.594 


3.343 


MRSA 


22.054 


6.651 


3.316 


CTEQ2M 


22.350 


6.708 


3.332 


CTEQ2MS 


21.662 


6.541 


3.312 


CTEQ2MF 


22.589 


6.788 


3.328 


CTEQ2ML 


23.357 


6.963 


3.354 


MRSD/— DIS 


22.190 


6.662 


3.331 


MRSH DIS 


22.404 


6.691 


3.348 


CTEQ2D 


22.670 


6.719 


3.374 



TABLE V. W and Z boson production cross section predictions, calculated using different pdf 
sets. All sets are in the MS, scheme except the last three, which use the DIS scheme. 



M w (GeV) 


a w (nb) 


a z (nb) 




80.05 


22.403 


6.671 


3.358 


80.23 


22.350 


6.708 


3.332 


80.41 


22.298 


6.745 


3.306 



TABLE VI. W and Z boson production cross section predictions, calculated for values of My/ 
one standard deviation below and above the world average (using the CTEQ2M pdf). 



The error on the mass of the W boson leads to an uncertainty in the value of the W 
cross section, and to a lesser extent the Z cross section (since sin 2 6 W is correlated with My/). 
The effect on the individual cross sections is small compared to that from the choice of pdf. 
However, in the ratio of cross sections, the two contributions are comparable. The effect of 
the My/ uncertainty is shown in Table [VJ. 

The last source of error considered for the variation of the calculated cross sections is 
the choice of factorization and renormalization scales. It is customary to set both scales 



equal to the same value f34|. We set the scales to the corresponding values of the vector 
boson masses. The uncertainty is estimated by varying the scales by a factor of two in either 
direction. The results are shown in Table [VlJ. The effect is small for the individual cross 
sections, as well as for the ratio. 

The effects of all the sources of error on the calculated cross sections are summarized 
in Table [VIII| . Using CTEQ2M for the central values, the theoretical predictions for the 
production cross sections at yfs =1.8 TeV are: 



a w = a(pp -v W + X) = 22 A + _\l nb, 
= a(pp ^Z + X) = 6.71±g;p nb, 
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Scale (GeV) 


a w (nb) 


a z (nb) 


o~w/o~z 


My/2 


22.259 


6.688 


3.328 


My 


22.350 


6.708 


3.332 


2M V 


22.421 


6.715 


3.339 



TABLE VII. Calculations of W and Z boson production cross sections for different values of 
the factorization and renormalization scales. My is the mass of the corresponding vector boson. 
(The CTEQ2M pdf and the nominal value of Mw are used.) 



Error Source 


Saw 


(nb) 


5a z 


(nb) 


5(a w /a z ) 


pdf choice 


+1.007, 


-0.688 


+0.255, 


-0.167 


+0.022, -0.020 


NLO pdf's 


+0.671, 


-0.671 


+0.201, 


-0.201 




5M W 


+0.053, 


-0.052 


+0.037, 


-0.037 


+0.026, -0.026 


Scale 


+0.071, 


-0.091 


+0.007, 


-0.020 


+0.007, -0.004 


Total Error 


+1.213, 


-0.967 


+0.327, 


-0.265 


+0.034, -0.033 



TABLE VIII. Summary of estimated uncertainties on the calculated W and Z boson production 
cross sections. The separate errors are added in quadrature to form the total error (assuming no 
correlation between error sources). 



— = 3.33 + 0.03. 
crz 

In order to compare the theoretical predictions of the W and Z boson production cross 
sections to experiment, it is necessary to multiply the cross sections by the vector boson 
branching fractions into the observed experimental channels. For the current study, only 
the electron and muon channels are used. 

Very precise values are available for the Z leptonic branching fractions. LEP measure- 
ments give p B[Z — > K) = (3.367 ± 0.006)%. For the W boson, we use a higher-order 
theoretical calculation |35j B(W —>■ Iv) = (10.84 + 0.02)%. Combining these branching frac- 
tions with the production cross sections quoted above gives the following predicted values 
for the cross section times branching fraction: 

a w B(W -> Iv) = 2.42t°;i? nb, 
a z B(Z -> ££) = 0.226l{j;£& nb. 



2. Results from Experiment 

The cross sections measured for W and Z boson production are calculated using the 
following formula: 
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FIG. 15. Measurements and predictions for W — > Iv and Z — > 11 cross sections. The results 
for this experiment are plotted as filled circles and those for the CDF experiment as open circles. 
The inner error bars represent the combined statistical and systematic uncertainties and the outer 
error bars include the uncertainty on integrated luminosity. The bands correspond to the range of 
predictions discussed in the text. 
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N 
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fbgdy 



Ae/Ldt 



where N b s is the number of events in our final data sample, fb g d is the fraction of the sample 
calculated to arise from background, A is the acceptance of the detector, e is the efficiency 
for accepted events to reach the final sample, and / Ldt is the integrated luminosity. 

The results are summarized in Table [TX[ Within the total errors, the measured cross 
sections are in good agreement with theoretical expectations. Our measurements are plotted, 
together with the predictions and other published experimental results 0] at ^Js = 1.8 TeV, 
in Fig. [L5. 



Channel 


W-ku/ 


Z e+e" 


W ^ \iv 


Z — ► 


N ohs 


10338 


775 


1665 


77 


Total Bgd(%) 


5.7 ±1.7 


4.0 ±1.4 


22.1 ± 1.9 


10.1 ±3.7 


Acceptance(%) 


46.0 ±0.6 


36.3 ±0.4 


24.8 ±0.7 


6.5 ±0.4 


£trig x e se i(%) 


70.4 ± 1.7 


73.6 ± 2.4 


21.9 ±2.2 


52.7 ±4.9 


f£dt (pb" 1 ) 


12.8 ±0.7 


12.8 ±0.7 


11.4 ±0.6 


11.4 ±0.6 


crB (nb)±(stat), 
(sys),(lum) 


2.36 ±0.02 
±0.08 ±0.13 


0.218 ±0.008 
±0.008 ± 0.012 


2.09 ±0.06 
±0.22 ±0.11 


0.178 ±0.022 
±0.021 ± 0.009 



TABLE IX. Observed cross section multiplied by the leptonic branching fraction for W and Z 



boson production. 
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C. Extraction of B(W -» tv) and T(W) from i? 

1. Phenomenological Considerations 

The leptonic branching fraction and the total decay width of the W boson can be ex- 
tracted from the measured ratio of the cross sections multiplied by the branching fractions 
of the W and Z bosons into leptons. The ratio R can be expressed as follows: 

_ a w B(W -> tv) _ a w 1 T(W -> tv) 

a z B(Z -> a z B(Z -> r(W) ' 1 ' ' 

Using an experimental result for R, the known B(Z — ► ££), and the prediction of <7w/°~z, 
a value for the leptonic branching fraction of the W boson follows: 

m r i») = w ; u> R = (001011 ±0 , 00011)fl , (7 . 4) 

Alternatively, using, in addition, a calculation of T(W — > ^i/), the full width of the PU 
boson width T(W) can be extracted: 



nW ) = { ° wl lf (W - (7.5) 



B(Z -> i2 
The leptonic width of the W boson can be written as: 



T(W -+lv) = — 7T (1 + ^ M ). (7.6) 
67TV2 



The corrections <5 SM have been calculated in the standard model by Rosner et al. [ 55| . Using 
G F = (1.16639 ±0.00002) x 10~ 5 GeV" 2 , M w = 80.23 ±0.18 GeV and 5 SM = -0.35% gives 
T(W — > £z/) = 0.2252 ± 0.0015 GeV, where the error is entirely due to the dependence on 
M w . 

In order to properly calculate the uncertainty on T(1U), it is necessary to take into 
account the correlation of errors on o~w/o~z and T(W — > tv) due to dependence on Mw- The 
product of these factors is shown in Table [X] for a one standard deviation variation in Mw- 
Taking the side with the larger variation as the error, the variation in the product is 0.0009 
GeV. The error on the product due to other sources is 0.0045 GeV; combining the errors 
in quadrature gives 0.0046 GeV. The product, using the nominal value of o~w/o~z = 3.33, is 
then 

—T(W -> tv) = 0.7499 ± 0.0046 GeV. (7.7) 
o~z 

Finally, using this value in the expression for T(W) leads to 

T(W) = (22.272 ± 0.137)-j-GeV. (7.8) 

R 
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M w 
(GeV) 


o~w/°~z 


T(W -> lv) 
(GeV) 


{a w /a z )T{W -> &/) 
(GeV) 


80.05 


3.358 


0.2237 


0.7512 (+0.0008) 


80.23 


3.332 


0.2252 


0.7504 


80.41 


3.306 


0.2267 


0.7495 (-0.0009) 



TABLE X. Calculation of the product (aw/o~z)F(W — > for a 1 a variation in Myy, using 
the pdf CTEQ2M. 



2. Result of Measurements 



The ratio of cross sections is given by: 

«(1 - O e z A z 



R 



Ng 8 (l-fl,)eW A w 



'bgdy 



(The dependence on the luminosity is completely canceled in the ratio.) Our results for e 
and fj, channels are: 

R e = 10.82 ± 0.41(stat) ± 0.35(syst), 
Rf, = 11.8ijj(stat) ± l.l(syst), 

and combined: 

R e+fl = 10.90 ± 0.52(stat©syst). 



This is consistent with previous measurements shown in Fig. [16 
Using this result, we obtain the branching fraction 

B(W -> lv) = (11.02 ± 0.52)%. 



(7.9) 



Combining this measurement with the calculation of the partial width of the W boson 
T(W — > lv), we obtain 



T{W) = 2.044 ± 0.097 GeV. 



(7.10) 



This is in excellent agreement with the prediction of the standard model, T(14 y ) = 2.077 ± 
0.014 GeV HH, and with the world average value, T(W) = 2.06 ± 0.06 GeV §. 

We can use our result to probe new possible decay modes of the W boson, such as 
decays into supersymmetric charginos and neutralinos |37j or heavy quarks [38]. Since our 
experimentally measured central value of r(W / )/r(H / — > lv) (the inverse of the branching 
fraction) falls below the mean predicted by the standard model, we use the asymmetric 
method to calculate limits on new decay modes M. From our data, we derive a 95% CL 
upper limit of 171 MeV on the width of unexpected decays of the W boson. If a new heavy 



quark exists, the limit for its mass is m q i > 61 GeV at the 95% CL (see Fig. [17]). Combining 
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UA1 (e+n) '83/'89 
UA2 (e) '88/'90 
Vs = 0.63 TeV <R> = 10.1 ± 0.6 



CDF (e) '88/'89 
CDF (n) '88/'89 
D0 (e) '92/'93 
D0 (n) '92/'93 
CDF (e) '92/'93 
Vs = 1.8 TeV <R> = 10.9 ± 0.3 



9 10 11 12 
R 



FIG. 16. Measurements of the ratio of the W — > Iv and Z — > 11 cross sections multiplied by 
their respective branching fractions. The results are shown as a function of the years of the data 
run. 



our result with other measurements [^S] gives a weighted average of = 2.062 ± 0.060 

GeV, and a 95% CL upper limit of 111 MeV on unexpected decays. 

Since the time that these results were first reported in a Letter || , the knowledge of the 
mass of the W boson has improved substantially. If we update the value used in Ref. |§ of 
M w = 80.23 ± 0.18 GeV to the current value of M w = 80.39 ± 0.06 GeV Q, the following 
results are obtained: 



a(W)B(W -> eu) = 2.35 ± 0.02 ± 0.08 ± 0.13nb, 



(7.11) 



B(W -> lu) = (11.03 ± 0.52)%. 



(7.12) 



T(W) = 2.054 ±0.097 GeV 
All other numbers reported change by much less than their uncertainties. 



(7.13) 



VIII. CONCLUSIONS 

D0 has measured the product of cross section and the lepton branching fraction for W 
and Z boson production in the electron and muon decay channels. We find 



42 




1 ' 6 20 40 60 80 100 

Heavy Quark Mass (GeV) 

FIG. 17. The width of the W boson as a function of a new quark mass. Our measurement is 
shown as a one standard deviation band with the central value represented by the solid line. The 
darker curve represents the prediction of the standard model as a function of quark mass. The 
short dashed line indicates the upper limit at 95% CL on the width of the W boson from our data. 

a(W)B(W -> eu) = 2.36 ± 0.02 ± 0.08 ± 0.13nb, 
a(W)B(W -> iiv) = 2.09 ± 0.06 ± 0.22 ± O.llnb, 
<r(Z)B(Z -> ee) = 0.218 ± 0.008 ± 0.008 ± 0.012nb, 

and 

<r(Z)B(Z -> tin) = 0.178 ± 0.022 ± 0.021 ± 0.009nb. 

Our values are in good agreement both with the 0(a 2 s ) QCD predictions using recent pdf 
sets, and with other measurements. 

Including theoretical calculations for crw/cz and B(Z — > 11), we measure 

B(W -»• \v) = (11.02 ± 0.52)%. (8.1) 

Adding the standard model prediction for T(W — > Zv), we find 

r(W) = 2.044 ± 0.097 GeV. (8.2) 

These results are in good agreement with the standard model, and allow us to set a limit 
on any new decay modes of the W boson. 
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